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ABSTRACT 



We identify the dual geometry of the hadron phase of dense nuclear matter and in- 
vestigate the confinement /deconfinement phase transition. We suggest that the low tem- 
perature phase of the RN black hole with the full backreaction of the bulk gauge field is 
described by the zero mass limit of the RN black hole with hard wall. We calculated the 
density dependence of critical temperature and found that the phase diagram closes. We 
also study the density dependence of the p meson mass. 
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1 Introduction 

Recently, related to the RHIC and LHC experiments, understanding strongly interact- 
ing QCD is requesting much attention. Although a powerful method for this subject, 
the lattice QCD, is being developed, when it comes to the dense matter problem, lat- 
tice calculation has difficulty and not much result is produced so far. The AdS/CFT 
correspondence [T] in the string theory, can shed many aspect of hadron theory [21 El HI O, 
El [TJ [HI [9l [To] as well as in strongly interacting quark gluon plasma [HI [121 [13 [HI [13 [IS] . 
The theory can easily accommodate the dense matter problem at least for deconfined 
phase [T71 [H [H [6] . However, for the the hadron phase, the status is not very clear since 
even the phase diagram is qualitatively different from that of the real QCD [H [18]. This 
may be traced back to the probe approach of the dense matter and one expect that if one 
fully account the back reaction of the gravity to the dense matter, one might overcome 
the situation. 

In fact, in the previous work[TO], it was suggested that if one considers bulk filling 
branes, one can easily account the full back reaction and the phase diagram actually 
closes. However, in that work, the back reacted geometry of the hadronic phase was not 
fully identified but was approximately treated as the thermal ads with electric potential. 

In this paper we consider the zero black hole mass limit of the charged black hole with 
hard wall installed as the low energy pair of the charged black hole. Although it has the 
naked singularity at the origin it is hidden in the wall and we do not find any physical 
difficulty. With this identification, density dependence of the physical quantities can be 
easily calculated. As a first example, we calculated how meson spectrum depends on the 
density of the baryonic medium. 
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The rest of paper follows: In the section 2, we will briefly review the dual geometry 
for the quark-gluon plasma and then propose what is the dual geometry for the hadronic 
phase. In the section 3, we will investigate the Hawking-Page transition. We will first 
study the fixed chemical potential case and then consider the fixed density problem. In 
the section 4, we will investigate the mass of the excited vector mesons in dense medium. 
In the section 5, we will summarize our results and discuss some future works. 



2 Dual geometry for QCD with quark matters 



In AdS/CFT correspondence, the boundary value of the bulk gauge field is coupled to the 
dual operator in the QCD side, which is the quark current. Furthermore the boundary 
value of the time-component gauge field is the chemical potential its dual operator is the 
quark number density. Our main interest here is to see how the critical temperature of 
the phase transition depends on the quark chemical potential. To describe the region for 
the high chemical potential, we need to consider the back reaction of the bulk gauge field 
in the dual geometry. Here, we will investigate the asymptotically AdS geometry dual to 
QCD with hadronic matters. 

We first review the gravity theory in the Mikowskian signature with introducing our 
conventions. The gravity action describing the five- dimensional asymptotic AdS space 
with the gauge field is given by 



Sm = J d^xV^ 



MN 



where 2k'^ is proportional to the five-dimensional Newton constant and is a five- 
dimensional gauge coupling constant. In the five dimensional AdS space, the cosmological 
constant is given by A = —G/R"^, where R is the radius of the AdS space. The equations 
of motion of this system becomes 

-,2 {^MpF^ - -^GmnFpqF^'^^ , 

= 9mv^G^^^G^«Fpq, (2) 
where M, TV = 0, 1, ■ ■ ■ , 4, x*^ = t and x'^ = z. Under the following ansatz 



n 



MN 



-GmnT^ + Gmn^ — — 
2 g- 



A, 
ds" 



^4 = (i = l,---,3), 
^ (-f{z)de + dxl + 



-dz' 



(3) 



the most general solution of Eq. 
(RNAdS BH) is 



known as the Reissner- Nordstrom AdS black hole 



1 — mz^ + (f'z^ , 



Ao = fi- Qz^ 



(4) 
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where the boundary space, denoted by = {x^jX*}, is located at 2; = 0. So, /i is a 
boundary value of Ao and Q is related to the black hole charge q through 



1' = ^Q'- 



(5) 



In the AdS/QCD context [TO], the gravitation constant and the five- dimensional 
coupling constant are related to the rank of the gauge group and the number of the 
fiavor Nf in QCD 

(6) 



2k2 87r2i?3 
so that Eq. ([5]) can be rewritten as 



Q 



and 



^2 Anm'' 



2Nf 



(7) 



2.1 Dual geometry of the quark-gluon plasma 

For investigating the Hawking-Page transition dual to the C/D phase transition in QCD, it 
is more convenient to consider the Euclidean version. Here, we summarize the Euclidean 
RNAdS BH shortly, which corresponds to the deconfinement phase described by the 
quark-gluon plasma. 

By the Wick rotation t —ir, the Euclidean version of the previous action in Eq. 
([1]) reads 



S = / d^xVG 



— {-n + 2A) + —^FmnF'''' 



where Gmn is the Euclidean metric ansatz 

i?2 



f{z)dT' + dx' + 



-dz' 



(9) 



The equations of motion for this system become 

K-2 

TZmn 



-GmnTI + Gmn^ 
2 g 



1 



FmpF^t — -GmnFpoF^'^ 



2 ^ - MPJ^N ~ -^^MN^ PQ± J 



= Smv^G^^^G^^Fpq. 
Under the following ansatz 

Ar 

A. 



(10) 



A{z), 

^4 = 0, 



(11) 



the most general solution is nothing but the Euclidean RNAdS BH, so that the metric 
factor f{z) is the same as one in the Minkowski version. Finally, the metric solution is 



ds' = ^ 



mz 



q'z'')dT' + d^' 



1 



1 — mz^ + q^z^ 



dz' 



(12) 
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From now on, we consider the Euclidean version only. Because F^t = —Ftz = dzAt, the 
Maxwell equation can be reduced to the simple form 

= dz(^^ dzA{z)^ , (13) 

and the solution is given by 

A{z) =i{fx-Qz^) . (14) 

Note that the imaginary number, i, in the above is very important to satisfy the Einstein 
equation in Eq. (fTOl) . which naturally appears due to the Wick rotation. 

From the metric in Eq. f[T^ . the outer horizon denoted by r+ should satisfy 

= f{z+) = 1 - mzl + 5^4- (15) 

Using the above, we can replace the black hole mass m with a function of the outer horizon 
2+ and the black hole charge q 

m=— + q^z\, (16) 
^+ 

which is useful for the later convenience. The Hawking temperature of the RNAdS BH is 
given by 

^™ -irS'- • 

For the norm of ||y4(z)|| = g'^'^A^Aj. at the black hole horizon to be regular, we should 
impose the Dirichlet boundary condition A{z^) = [lOl HI [5l [191 [20] , which gives a relation 
between Q and /i 

Q' = ^. (18) 

4 

Inserting this relation and Eq. ([5]) into Eq. f[T7j) . we can find clS db function of /i and 

Trn 

= w U + ^ ~ "^^^ j • ^''^ 

To describe a system having the fixed chemical potential, we should impose the Dirich- 
let boundary condition, A{0) = ifi, at the UV cut-off z = e where e is very small. Then, 
the on-shell action becomes 

where V3 corresponds to the spatial volume of the boundary space. In the above action, 
the superscript D and the subscript RN imply the Dirichlet boundary condition at the 
UV cut-off and the RNAdS BH, respectively. Since this action has a divergent term 
when e ^ 0, we should renormalize it. For the renomalization, we use the background 
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subtraction method, in which the on-shell action for the AdS space is subtracted from 
S^jtf. The on-shell action for the AdS space is 



SAds = ^ j dT r dz VG(-n + 2K) 



After identifying the circumference of the RNAdS BH and the AdS space at the boundary, 
we can rewrite /3 as 

/3 = ^ (l - \me' + 0(6«)) . (22) 
Therefore, the on-shell action of the AdS space becomes 

^^'^^-^TWV^ 24 3^4^' ^^^^ 

where Eq. ([5]), Eq. (fT6|) and Eq. (ITSl) are used. As a result, the renormalized action of 
the RNAdS BH is given by 



-l-RN 

and the grand potential becomes 



2z\ Sg'^R^z^^ 



(24) 



RN — Orn-'-RN 

'3/1 ^2 ,,2 



k2 \2zl Sg^R^z^^ 

For describing the dependence of the particle number in this system, we should consider 
the free energy F obtained by the Legendre transformation of the grand potential. In the 
thermodynamics, it is given by 

F = + /iiV, (26) 

where N = Using the fact that 2;+ can be represented as a function of /i and T like 

Eq. (HM . N is given by 

N = -iQVs, (27) 

where the particle number is proportional to the quark number density Q. As a result, 
the free energy becomes 

9 n 

F = n + —fiQVs. (28) 
9^ 

Interestingly, this result can be reobtained from the action Eq. ([8]) with the different 
boundary condition at the UV cut-off. For changing the Dirichlet boundary condition 
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into the Neumann boundary condition, we should add a boundary term to fix Q. Then, 
the renormahzed action is given by 



^RN — SrN + ^bi (29) 

where the superscript imphes the Neumann boundary condition at the UV cut-off and 
the boundary action Sh is given by 

Sb = \ [ d^xVm n^'A^^FMN. (30) 

9 JdM 

In the above, the unit normal vector is given by n^'^ = |o, 0,0,0, ^\'^T(^^ and G*^^^ = 

^f{z) is a determinant of the four- dimensional boundary metric. Using the solution for 
the bulk gauge field in Eq. ( |T4l) . the boundary term becomes 



So, the free energy reads 



Sb = ^^I^Q. (31) 
Trn 9 



- 9 R 

F = S^j^TjiN = ^ + ^I^QVs, (32) 



which is the same as one obtained from the thermodynamics in Eq. (l28il . From these 
results, we may conclude that the five- dimensional bulk action with the Dirichlet or 
Neumann boundary condition at the UV-cut off corresponds to the grand potential or 
free energy of the dual QCD, respectively. 



2.2 Dual geometry of the hadronic phase 

In QCD, it is well known that there exist the hadronic or confinement phase at the low 
temperature. In the absence of quark matters, the Schwarzschild AdS black hole (SAdS 
BH) corresponds to the deconfinement phase. At the low temperature, the dual geometry 
for the confinement phase is given by the thermal AdS with the IR cut-off, which is 
needed to explain the confining behavior [21]. What is the dual geometry describing the 
hadronic phase, in other words the confinement phase with quark matters? As shown in 
the previous section, we should include the bulk gauge field to explain the quark matters. 
Therefore, the geometry corresponding to the hadronic phase has to be a deformed AdS 
including the backreaction of the bulk gauge field, which is not a black hole. The metric 
we find out to answer the above question, is 

ds' = ^ f (1 + qh')dT' + dx' + -^—^dzA , (33) 
\ 1 + q^z*^ J 

which together with Eq. (HM satisfies the Einstein and Maxwell equation and becomes 
a AdS space asymptotically. This solution can be also easily obtained from the RNAdS 
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BH by taking m = 0. For the convenience, we call this solution as a thermal charged 
AdS (tcAdS) solution. Especially, in the case of g = 0, the tcAdS and RNAdS BH are 
reduced to the thermal AdS (tAdS) and SAdS BH, respectively. Here, our proposition 
is that the tcAdS having the IR cut-off zm is the dual geometry corresponding to the 
hadronic phase. 

Now, we study the thermodynamics of the tcAdS with the fixed chemical potential, 
for which the Dirichlet boundary condition is needed at the UV cut-off. At the IR cut-off 
we impose another Dirichlet boundary condition 

A{zir) = iafi, (34) 

where a is a arbitrary constant and will be determined later. This IR boundary condition 
together with Eq. (fT4|) gives a relation between fi and Q 

Q = <1^. (35) 

^IR 

Using this, the on-shell action of the tcAdS is given by 

^ VsR^ 1/1 1 2k^ (l-a)V 



- T,. [e^ ?,gm- zj^ ) ' ^^^^ 

where the subscript tc means the tcAdS. To renormalize this action, we subtract the 
on-shell action of the AdS space in Eq. (12 ip . with the identification between the circum- 
ferences of two backgrounds at the UV cut-off 

= (1 + 0{e')) . (37) 

J-tc 

Then, in the limit of e — the renormalized action of the tcAdS becomes 
.D _ V,R' 1 f 1 2k' (1-«)V \ 



k' Tt, \zj^ 3gm' zjj, J 

From this, the particle number is given by 

9 op 

N = -il ~ a)^QV^. (39) 
3 g' 

Since the action with the Neumann condition should be a free energy as previously men- 
tioned, iiN = SbTtc should be satisfied. Using the fact that the boundary action Sb of the 
tcAdS is given by 

a should becomes —1/2 for the consistency. As a result, the renormalized action becomes 

V^R^ I f I , 3^2 



Sb=^—QVs, (40) 



with the following relation 

1^ = \qAr- (42) 
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3 Confinement /deconfinement phase transition 



In QCD, there exists the C/D phase transition, which is dual to the Hawking-Page transi- 
tion in the gravity theory side. So it is an interesting question to ask how the C/D phase 
transition depends on the quark matters. As mentioned previously, we should introduce 
the IR cut-off to describe the confining behavior of the hadronic phase, which is called 
the hard wall model. There were many interesting works to explain QCD depending on 
quark matters. In our previous work |9] , we studied the dependence of the quark (num- 
ber) density in the C/D phase transition by considering the gauge field fluctuation on the 
tAdS and SAdS BH backgroud. Unfortunately, its result is valid only in the regime of the 
low quark number density so that it can not explain the dependence of the high chemical 
potential and quark density. Anyway, the crucial point of the tcAdS and RNAdS BH 
backgrounds is that since they include the full backreaction of the gauge field, we can 
investigate the dependence of quark matters even in the high chemical potential or quark 
density regime using this model. 



3.1 Fixed quark chemical potential 

For the fixed chemical potential, to describe the Hawking-Page transition we calculate the 
difference between the on-shell actions of two backgrounds, with the Dirichlet boundary 
condition at the UV cut-off, 

AS = S^^ - Sg, (43) 

with 

qD 

'-'RN - 





1 


(--- 


2k^ /i2 




Trn 


[e^ 4 




VsR' 




'1 1 




k2 


Ttc ( 


.~'~Wr 


2^2/22 ,2^ 



k2 T,Se' 2gm'zj^)- ^^^^ 

After requiring the same circumference of r at the UV cut-off, the difference becomes 

. . ^ VsR^J_ (I L : J^iLZ _ r4^^ 

TnN V4« 24 ^ 2gm^ z]^ -ig^^ zl) ' ^ 

which is the same as one obtained from the renormalized actions in the limit of e ^ 0. 
The Hawking-Page transition corresponding the C/D phase transition occurs at AS* = 0. 
Note that the C/D phase transition occurs only in the range of 2+ < zir. So we will 
consider the case of z^^ > zir from now on. Suppose that AS* is zero at a critical point 
Zj^ = Zc- In the case of z+ < Zc, AS becomes negative. So the RNAdS BH is stable, 
which implies that the dual boundary theory is described by quark-gluon plasma or the 
deconfinement phase. In the Schwarzschild black hole case, since there is no black hole 
charge, the dual theory is described by the gluon only, without including quark matters. 
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Figure 1: The deconfinement temperature depending on the chemical potential, which does not 
contain the quark mass. 



Anyway, for Zc < Zj^ < Zm the stable space is the tcAdS. Since the tcAdS corresponds to 
the confinement phase, the dual QCD describes the hadronic matters. 

For the later convenience, we introduce dimensionless variables as the following 

~ = ^ 

ZiR 

j^c = f^cZlR, 



TcZiR, 



(46) 



where the subscript c means the critical values representing the C/D phase transition 
point. Then, the critical chemical potential flc and temperature Tc can be represented as 
functions of Zc 




2~zt) 



zl{9~zl-2) 



Tz 



9z 



~2 



(47) 



To obtain a well-defined chemical potential, the inside of the square root in the above 
should become positive. So the allowed range of Zc is give by 0.4714 < 5^ < 0.8409. 
Though at flc = the deconfinement temperature Tc = 122 MeV is lower than the lattice 
estimation = 175 ~ 190 MeV [221 123] as noted by Herzog [21], we use = 122 MeV 
and look at the qualitative behavior of the deconfinement temperature depending on the 
chemical potential. We obtain the expected phase diagram, which is given in the Figure 
1. Notice that the deconfinement temperature decreases more quickly as the ratio Nf/Nc 
goes up. In the holographic models with probe brane approach, the phase diagram does 
not close. However, in our case the phase diagram closes due to the backreaction of the 
gauge field as emphasized in Ref. (TU]. In the zero temperature the critical value of the 
phase transition, fic(Tc = 0) is given by 1002 MeV for Nf/Nc = 1- Note that our definition 
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Figure 2: The Hawking-Page transition temperature depending on the quark number density. 

of the chemical potential does not include quark mass. From this, we can evaluate the 
the critical baryon chemical potential with mass, 

fiB = 3iJ + TUB ^ 4:GeV , (48) 

where is the baryon mass. For more realistic value, we should use different normal- 
ization or other models like the one suggested in Ref. |24j . 



3.2 Fixed quark number density 



For the fixed quark number density, we should add a boundary term for fixing Q, which 
corresponds to the Neumann boundary condition at the UV cut-off. The on-shell actions 
of the RNAdS BH and the tcAdS are given by 



S 



N 
RN 



N 
tc 



Trn 



-tc 



1 



1 

-'IR 



(49) 



Then, the difference between two on-shell actions becomes 



qN 

'-'rn 



RN 



''IR 



2z, 



(50) 



Like the previous section, we assume that Zc is a critical value of z+ where the C/D phase 
transition occurs. For z^ < Zc, the RNAdS BH is more stable than the tcAdS. On the 
contrary, the tcAdS is stable for 2+ > Zc- 

After introducing a new dimensionless variable, Qc = QcZ^^.^ the dimensionless quark 
number density and critical temperature are given as functions of Zc 



Qc 



{l-2zt) 
2Nf~zt{h~zl-2) 
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Figure 3: The p meson mass depending on the chemical potential. 



1 



\ _ ~4 (1 - 'izt) 



2 (552 -2) J 



(51) 



5 < < which is 



In the first relation, Qc is well defined only in the range, 
0.6325 ^ 5c < 0.8409. Using Eq. (15T|) . we numerically draw the deconfinement tem- 
perature depending on the quark density in the Figure 2. As shown in the Figure 2, 
the deconfinement temperature decreases as the quark density increases. Furthermore, 
when the number of fiavor becomes large, the deconfinement temperature decreases more 
quickly than one having the smaller number of fiavor. 



4 The mass of the excited vector mesons 

At first, we consider the p meson mass at zero temperature and finite baryon chemical 
potential p,Bi which is given by three times as much as the quark chemical potential, 
/i^ = 3/iq. To describe this, we should consider the fiuctuation of the vector field in the 
tcAdS. The ansatz for the vector field fluctuation is given by 

5A^ = V,{z,p)e'P-\ (52) 

Since the rotation symmetry in the r-Xi plane is broken, we should distinguish the 5Aq 
with 5Ai. The equation of motion for 5Aq is the exactly same as one obtained in tAdS. 
So, from now on we will concentrate on the 5Ai fluctuation having the following equation 
of motion 

= dlV. - - \ \ ^ d,V, + mlV,, (53) 

where the meson mass is given by = —p"^. After introducing the new dimensionless 
variables 

~ z 

Vi = ViZiR , z = and fhm = rrimZiR, (54) 

^IR 
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Figure 4: For Nf/Nc = 2/3, the chemical potential dependence of the excited meson masses. 

Eq. fl53l) can be rewritten as 

= dIV, - d,V, + (55) 

For g = 0, the above equation reproduces the result in tAdS [2Tj. Due to the C/D phase 
transition, the range of the quark chemical potential fx at the zero temperature is limited 
from to 5.373 for -jf = \- Similarly, jl runs from to 3.799 for ^ = § and from to 
3.10217 for Nf/Nc = 1- Through the numerical evaluation, we find the relation between 
the chemical potential and the p meson mass in the Figure 3. As shown in the figure, 
the p meson mass decreases as the quark chemical potential increases. Note that in the 
Figure 3, the right ends of the curves implies the C/D phase transition point. In addition, 
the larger the number of flavor becomes, the more quickly the p meson mass decreases. 





/i = 


p = 0.245 


p = 0.491 


p = 0.736 


p = 0.982 


p = 1.227 


mass of the 1st 


0.774 


0.724 


0.622 


0.530 


0.458 


0.404 


mass of the 2nd 


1.775 


1.737 


1.702 


1.704 


1.724 


1.750 


mass of the 3rd 


2.782 


2.758 


2.743 


2.747 


2.755 


2.762 



Table 1: The meson masses of the excited states depending on the chemical potential in 
the GeV unit. 



Finally, we further calculate the masses of the higher excited meson states. After the 
numerical evaluation, we draw several excited meson masses for ^ = | in the Figure 4 
(see also Table 1. for the numerical values of the masses in the several points of p). For 
the p meson, its mass decreases as the chemical potential increases. In the higher mode 
cases, their masses decrease for some period of p and then increase as p increases (see 
Table 1). 
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5 Discussion 



In this paper, we proposed the thermal charged AdS space as the dual geometry corre- 
sponding on the hadronic phase of QCD, which is the zero mass limit of the Reissner- 
Nordstrom AdS black hole with the hard wall. This tcAdS was installed as the low 
temperature pair of the RNAdS black hole. By comparing the on-shell actions of two 
backgrounds, we investigated the confinement /deconfinement phase transition depending 
on the chemical potential or the quark number density. Interestingly, we found out that 
above the critical chemical potential there exists only the deconfinement phase even at 
zero temperature and evaluated the critical chemical potential depending on the flavor 
number. In addition, using the dual geometry of the hadronic phase, we calculated how 
the mass spectrum of the vector mesons depends on the chemical potential in the baryonic 
medium. 

In this paper, we evaluated the baryonic chemical potential (~ 4 GeV), which is 
too high comparing with the known result 1 Gev). So it is important to know how 
to cure this discrepancy. Another interesting problem is to study the chemical potential 
dependence (or the quark density dependence) of other physical quantities. We will report 
those results elsewhere. 
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